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ABSTRACT 

0^ ■ The spectra of the first galaxies and quasars in the Universe should be strongly 



0^ 



absorbed shortward of their rest-frame Lya wavelength by neutral hydrogen (HI) in 
the intervening intergalactic medium. However, the Lya line photons emitted by these 
C I sources are not eliminated but rather scatter until they redshift out of resonance and 

"^ ' escape due to the Hubble expansion of the surrounding intergalactic HI. We calculate 

the resulting brightness distribution and the spectral shape of the diffuse Lya line 
emission around high redshift sources, before the intergalactic medium was reionized. 
^ . Typically, the Lya photons emitted by a source at Zg ~ 10 scatter over a characteristic 

^^ , angular radius of ~ 15" around the source and compose a line which is broadened 

and redshifted by ~ 10'^ km s^^ relative to the source. The scattered photons are 
CN I highly polarized. Detection of the diffuse Lya halos around high redshift sources would 

Qs I provide a unique tool for probing the neutral intergalactic medium before the epoch of 

On ■ reionization. On sufficiently large scales where the Hubble fiow is smooth and the gas 

is neutral, the Lya brightness distribution can be used to determine the cosmological 
mass densities of baryons and matter. 
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^ ■ Subject headings: cosmology: theory - line: profiles 
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1. Introduction 

Following cosmological recombination at z ~ 10^, the Universe became predominantly neutral 
and hence optically-thick to Lya photons. The first galaxies that lit up at lower redshifts were 
surrounded by a neutral intergalactic medium. Their observed spectrum should therefore show a 
deep trough shortward of their rest-frame Lya wavelength due to absorption by neutral hydrogen 
(HI) along the line-of-sight, the so-called Gunn- Peterson effect (Gunn & Peterson 1965). For 
typical cosmological parameters, the optical depth at the Lya resonance is so large that one might 
naively expect the damping wing of the Lya trough to eliminate any trace of a Lya emission line in 
the observed source spectrum (Miralda-Escude & Rees 1998; Miralda-Escude 1998). However, here 
we point out that the Lya line photons absorbed by intergalactic HI are subsequently re-emitted 
and hence do not get destroyed^. Rather, these photons scatter and diffuse in frequency to the red 



^At the redshifts of interest, Zs ~ 10, the low densities and lack of chemical enrichment of the IGM make the 
destruction of Lya photons by two- photon decay or dust absorption unimportant. 



of the Lya resonance due to the Hubble expansion of the surrounding HI. Eventually, when their 
net frequency redshift is sufficiently large, they escape and travel freely towards the observer. In 
this paper we calculate the resulting brightness distribution and spectral line profile of the diffuse 
Lya radiation around high-redshift sources^. 

The lack of a Gunn- Peterson trough (i.e., the detection of transmitted flux shortward of the 
Lya wavelength at the source redshift) in the observed spectra of galaxies at z ^ 5.6 (Dey et al. 
1998; Hu, Cowie, & McMahon 1998; Spinrad et al. 1998; Weymann et al. 1998) imphes that 
most of the intergalactic hydrogen in the Universe was reionized before then. Indeed, popular 
cosmological models predict that reionization took place around a redshift 2; ~ 10 (e.g., Haiman &: 
Loeb 1998a,b; Gnedin & Ostriker 1997). At earlier times, the abundance of ionizing sources was 
small and each of these sources produced an expanding HII region in the surrounding intergalactic 
medium (IGM). The volume filling factor of HII increased dramatically as the number of ionizing 
sources grew and their associated HII regions expanded. Eventually, these HII regions overlapped 
and hydrogen throughout most of the cosmic volume became ionized^ (Loeb 1997). Subsequently, 
the ionizing background radiation penetrated the denser HI condensations around collapsed 
objects (Miralda-Escude et al. 1998; Barkana & Loeb 1999). 

The physics of reionization involves complicated gas dynamics and radiative transfer. 
Numerical simulations had only recently started to incorporate the relevant equations of radiative 
transfer rigorously (Abel, Norman, &: Madau 1998). The detection of diffuse Lya photons around 
high-redshift sources offers an important empirical tool for probing the neutral IGM before and 
during the reionization epoch, and can be used to test related theoretical calculations. In this 
context, one is using the Lya source as a light bulb which illuminates the HI fog in the surrounding 
IGM. Fortunately, the highest redshift galaxies are found to be strong Lya emitters (see, e.g. Dey 
et al. 1998), probably due to their low dust content. 

In this paper we calculate the intensity distribution of the diffuse Lya line in the simplest 
setting of a pure Hubble flow in a neutral IGM around a steady point source. We leave more 
complicated configurations for future work. In §2 we describe the formalism and derive analytical 
solutions in the diffusion regime. In §3 we present a complete numerical solution to this problem 
using a Monte-Carlo simulation. The detectability of Lya halos is discussed in §4. Finally, §5 
summarizes the implications of our results. 



^The photons absorbed in the Gunn- Peterson trough are also re-emitted by the IGM around the source. However, 
since these photons originate on the blue side of the Lya resonance, they travel a longer distance from the source than 
the Lya line photons do, before they escape to the observer. The Gunn-Peterson photons are therefore scattered from 
a larger and hence dimmer halo around the source. The Gunn-Peterson halo is made even dimmer relative to the 
Lya line halo, by the fact that the luminosity of the source per unit frequency is often much lower in the continuum 
than in the Lya line. We therefore focus on the Lya line halo in this paper. 

■^Note, however, that even after much of the IGM was ionized, the optical depth at the Lya resonance was still 
substantial due to the small residual abundance of HI. Estimates indicate that the reionized Universe provided a 
measurable transmission of Lya photons only around z <.7 (Haiman & Loeb 1998; Miralda-Escude et al. 1998). 



2. Analytical Formalism 

Although the radiative transfer of line radiation was treated in the past for different geometries 
of stationary (Harrington 1973) or moving (Neufeld & McKee 1988) atmospheres, no particular 
attention was given to the cosmological context. As a first treatment of this involved problem, we 
consider here the transfer of the Lya line in a uniform, fully neutral IGM, which follows a pure 
Hubble expansion around a steady point source. As shown later, the diffuse Lya emission around 
a high-redshift source could extend over sufficiently large radii where these idealized assumptions 
are indeed valid. 

For newly created Lya photons the IGM is very opaque, and during scattering events the 
photon frequencies are redistributed symmetrically away from the line center due to the isotropic 
distribution of thermal velocities of the hydrogen atoms. However, as the photon frequencies drift 
away from resonance and the medium becomes less opaque, the asym,m,etric redshift bias imposed 
by the Hubble expansion of the surrounding IGM becomes dominant. Since the IGM is still highly 
opaque when the redshift effect starts to dominated, we focus in the following only on the Hubble 
expansion when calculating the observed intensity distribution around the source. In terms of 
the comoving frequencies used here, this implies that the scattering may be regarded as coherent 
(elastic). 

We consider a source surrounded by a spherically-symmetric atmosphere of neutral hydrogen 
with a linear radial velocity: v{r) = Hgr, where Hg is the Hubble expansion rate at the source 
redshift, Zg. Let / = I{v,ii,r) be the specific intensity (in photons cm^^ s^^ sr^^ Hz^^) at 
comoving frequency u m. & direction /i = cos 9 relative to the radius vector at radius r, as seen by 
an observer comoving with the gas (i.e. in the cosmic rest frame). Assuming isotropic coherent 
scattering, the comoving transfer equation for a line at a resonant frequency vq is then given by, 

/u— + ■— + a—=Xu{J-I) + S, 1 

or r Ofi Oh' 

where ly is the frequency redshift, namely the resonant frequency vq minus the photon frequency; 
Xi^ is the scattering opacity at v; J = ^ J_^ dfi I is the mean intensity; S is the emission function 
for newly created photons (in photons cm~^ s~^ sr~^ Hz~^); and 



a ^ :^ = 2.7 X W-'-'ho 



Om(1 + z^f + (1 - Qm - f^A)(l + Zsf + Qa] ^''^ Hz cm"\ (2) 



with Q^M and f^A being the density parameters of matter and the cosmological constant, and h^ 
being the current Hubble constant, iJoi in units of 100 km s~^ Mpc~^. Note that at sufficiently 
high redshifts the value of a does not depend on Qa- The source function on the right-hand-side 



*The neutral IGM is even colder than the microwave background prior to reionization. The corresponding thermal 
velocities of hydrogen atoms at jz ^ 30 are smaller by ^ 3 orders of magnitude than the Doppler velocity shift required 
for the cosmological escape of resonant Lya photons {^^ 10'^ km s~^). 



of equation (|lj) can be written as 



Na 



where Na = const is the steady emission rate of Lya photons by the source (in photons s~^). 
In terms of our frequency variable, the Lya opacity is given by (Peebles 1993, p. 573), 

V Z/Q/ i^ + A^[l — z^/z^oj /lOTT^ 

where vq = 2.47 x 10^^ Hz is the Lya frequency, A = 6.25 x 10® s^^ is the rate of spontaneous 
decay from the 2p to the Is energy levels of hydrogen, and 

/ 3r2A2 \ 
/3 = ( ^^^^3^ J nm = 1.5 QbhUl + z,f cm-^ Rz\ (5) 

Here nni is the mean hydrogen density in the (neutral) IGM at a redshift z^, expressed in terms of 
the current baryonic density parameter, Jl^, and the normalized Hubble constant, /iq. For typical 
cosmological parameters, the mean-free-path at the line center is negligible compared to the size 
of the system. Hence, a distant observer would see only those photons which scatter to the wings 
of the line, where to a good approximation the opacity scales as, 

Xu ~ ^, (6) 

with (3 being a constant due to the assumed uniformity of the IGM. This approximate relation 
holds as long as the frequency shift is small, z^/z^o ^ 1) and is adequate throughout the discussion 
that follows. 

Before proceeding, we should point out that the transfer problem described by equation (||) 
is completely analogous to that of time-dependent monochromatic isotropic scattering, where the 
frequency v here plays the role of "time." This is easily understood, since a photon is subjected to 
a constant redshifting due to the Hubble expansion, and the scattering is coherent, so its frequency 
is a perfect surrogate for time. This implies that the problem can be viewed as an initial value 
problem in frequency, so the solution for a fixed frequency determines the behavior of the solution 
at all other frequencies. One should note, however, when using this analogy that the opacity law 
here is a decreasing function of "time." 

We would now like to normalize the frequency shift and radius in the problem by convenient 
scales. An appropriate frequency scale, z/^, is introduced by equating the Lya optical depth from 
the source to the observer to unity, 

h {Vi, + ar) OLV^ 

yielding 

z/^ = - = 5.6 X lO^^f^fc/io ff^Af (1 + 2s)"^ + (1 - ^M - ^a)(1 + 2s)"^ + ^a(1 + zs)"1 '^''^ Hz. (8) 
a L J 



Note that for the popular values of Q^ ~ 0.05, JIm = 0.4, JIa = 0.6 and /iq = 0.65 (e.g., 
Turner 1999, and references therein), the frequency shift of the photons which escape to 
infinity is (u^/uq) ~ 4 x 10~^ at Zg = 10. This corresponds to a spectroscopic velocity shift of 
{y-i,/vQ)c ~ 10^ km s~^ at the source. The proper radius where the Doppler shift due to the Hubble 
expansion produces the above frequency shift is given by r^ = z^^/q, or 

/? 2.1 X 10^5 (Qb/^M) cm ( Vt 



o? 



1 + (1 - fij, - !iA)SJi/(l + z,)-i + (nA/fij;)(l + z,)-'' 



"•^n;^)"-"- <^> 



The last equality was obtained for high redshifts, z^^ \. We thus find that for an ^m = 1) ^A = 
cosmology or at high redshifts for any other cosmology, the physical distance r^ is independent of 
the source redshift and depends only on the baryonic mass fraction, Fb = Qb/i^M- This fraction 
can be calibrated empirically from X-ray data on clusters of galaxies (Carlberg et al. 1998; Ettori 
& Fabian 1999) which yield a value of Ff, ~ 0.1 for h^ = 0.65. Note that v^, was derived at the 
source rest frame; its value is lower by a factor of (1 + Zg) in the observer's frame. 

It is convenient for further developments to rescale the variables in the problem and make 
them dimensionless. Radius and frequency are rescaled using the characteristic quantities r^ and 
z^j,, that is, f = (r/r^) and P = [v/v^). We rescale all radiation quantities using the characteristic 
quantity I^, = Na/i^i^v*), that is, / = I /li,, J = J/Ii,, etc. We also define the rescaled source 

term, 

Sr^^ 1 

h (47r)2f 

Equation (|l|) can then be written in dimensionless coordinates. 



s = ^ = Tj^^2^mm (10) 



di_ (^^^af ^_}_(j_t\ q 

df f dfj, du u"^ ^ ' 

Next, we take the first two angular moments of equation (|ll|). In addition to J, these 
will involve the two moments K = -^ J_^ dfi fi^I, and H = ^ /-i dfi fiL In order to close the 
above system of equations, one needs a relation between J and K. The closure relation is often 
parameterized through the Eddington factor / = /(i/, r) = K/J (Hummer & Rybicki 1971; 
Mihalas 1978). This parametrization yield for the comoving transfer equation (Mihalas, Kunasz, 

& Hummer 1976), 

1 f)(r'2 m fIJ 

(12) 



1 dirm) 1 5J _ ^ 

^2 Qf Qy 




d{fJ) 1 {'if -I) J 1 dH _ 
df f " ' di> 


H 



and 

fi(fh (^f - ^^ ^ riff M 

(13) 



2.1. The Extended Eddington Approximation 

We can close the above set of equations by specifying the radial and frequency dependence 
of the Eddington factor /. The simplest choice is to use the Eddington approximation / = 1/3, 
which is expected to be valid at large optical depths in the medium. However, for the moment we 
assume / to be a constant, but not necessarily equal to 1/3. 

For convenience, we define the rescaled variables h = r^H, j = f'^J, and s = r'^S. This yields, 

|^ = -|i + ». (14) 

or ou 



fit-^-m = -i^.-To- (15) 



and 

or r ov v^ 

By taking the f~derivative of equation ( p^ and substituting the i>~derivative of equation (|14[) into 
it, we get 



.. ,/7p-[l-/]^ =Tp^ + ^Tp-Tp-^. (16) 

or \ Or r J ow'' v^ ov ov v^ 

If we now define the function 

g = ff-[l~f]/fj ^ ffW-n/fj, (17) 

then equation (|lq) can be written as, 

^ _ fp2 ( 1 A,~[l-/]//gg _ 1^^ = p2^ ~ .18^ 

dv ^ \f[^-f]/fd/ dr fdDy 5i> ' ^ ' 

where s = /f PZ-i]//^. 

Equation ([iq ) resembles the causal diffusion equation (Narayan, Loeb, Sz Kumar 1994, and 
references therein), with the substitution of frequency shift for time, and a "diffusion" coefficient 
D = fiP'. The Lya photons are expected to redshift away from the line center at a rate which is 
dictated by this diffusion coefficient. Inspection of the second-order terms provides the condition 
for "causality" in this problem, namely v > fjyfj. This condition has a simple interpretation: 
the frequency shift at any given radius is greater than the Doppler shift due to the Hubble 
velocity there, since the photon could have scattered multiple times before arriving at that radius. 
Unfortunately, the line opacity is not constant and the above diffusion equation is complicated by 
the frequency dependence of the diffusion coefficient. 



Equation ( p!8| ) requires a particular choice for the value of /. There are two physical regimes 
where the assumption of / = const holds: (i) the diffusion regime where the optical depth is high 
and / ~ 1/3, and (ii) the free-streaming regime where the optical depth is low and / ~ 1. Note 
that g = 3-/ in the first regime and g = f^J in the second. In the diffusion regime, the number 
of scatterings is large and so v ^ ar or v ^ f. This implies that the d^g term can be neglected 
relative to the V^g term in equation (18) and the problem is simplified considerably. Such a 



simplification is not possible in the free-streaming case. We therefore focus next on the diffusion 
regime where f <^ v <^ 1, and / ~ 1/3. 



2.2. The Diffusion Solution 



In the diffusion regime, equation (18) reads 

y'"^ = mT£^,- (19) 



dJ D^ 1 d odJ ,, , Sir) 

r — = div) —^ . 

du 3 f2 df df (47r)2f2 



With the definition of a new variable, a = (P'^/9), we get the standard diffusion equation in 

spherical geometry, 

f) 1 1 

--VV=-<^(.)<^(f), (20) 

where r is the radius vector, and where a now plays the role of the usual time variable. The 
solution to this equation for u > Q (redshifted photons) is 

This solution satisfies the integral relation, 

{Anf / Jr'^df = 1, (22) 

Jo 



which follows by performing the Gaussian integral, or directly from the diffusion equation (|19D. In 
a way this may be viewed as a result of photon conservation, since frequency here is analogous to 
time in the usual diffusion equation. 

Eventually the diffusion solution becomes invalid when the opacity is small enough to allow 
photons to escape to the observer. Most of the Lya photons escape to infinity from a smaller 
radius than r* as they scatter many times and hence redshift in frequency more than ar at any 
given radius r. Even though z/j, and r^ may not represent the true parameters of the escaping 
photons, equations (0) and (pi) provide the scaling of the exact numerical solution with model 
parameters, since we have already shown that these constants can be used to normalize the basic 
transfer equation (|I]) and bring it to a scale-free form. Next we proceed to derive the full numerical 
solution for the escaping radiation in the dimensionless variables u and f. 



3. Results from a Monte-Carlo Simulation 

The lack of an analytical solution applicable to both the diffusion and free-streaming regimes 
motivated us to find solutions by numerical means. The Monte Carlo method suggested itself 
by its generality and ease of formulation. The usual disadvantage of this method of having to 
treat large numbers of photons is largely irrelevant in the present case, because all the physical 
parameters have been scaled out of the equations, and only one Monte Carlo run is needed to 
solve the problem completely. 



The Monte Carlo method follows the histories of a large number of photons in accordance 
with the following rules: 

Step 1: A newly created photon is chosen from the distribution given by the form of the emission 
function S. Ideally we would like this to be a photon with f = and i> = 0. However, because of 
the singularity of the opacity, this cannot be done, so in practice we choose our new photons with 
a small, but nonzero, frequency (e.g., D = 10~^) and use the diffusion solution (E^) to determine a 
random starting value of f. 

Step 2: Suppose a photon has been scattered (or has been newly created) at a radius f and with 
a frequency D. A random direction is then chosen from an isotropic distribution; this is most easily 
done by taking /i = 2R — 1 where i? is a random deviate uniform on the interval (0, 1). 

To choose a step length for the photon flight, we use the fact that the optical depth r along 
the ray to the next scattering follows the exponential distribution law, exp(— r). This can be 
simulated by chosing r = —InR', where R' is another random deviate uniform on the interval 
(0, 1). To convert this to distance, we compute r as, 

^ r^ dl' 1 1 

X,d£'= , =£-^-,, (23) 

io (P + f)2 i^ D + i ^ ^ 

where i is the pathlength along the ray, and i = i/r-^. Taking i -^ cxd, we see that there is a 
maximum optical depth to infinity, 

''"max — ~- \^^) 

V 

If r is greater than Tmaxj the photon escapes and we proceed to Step 3. If r is less than Tmax; a 
scattering event will occur. We may then solve equation (|23|) for the pathlength. 



1= ""_ . (25) 

''"max ''" 

At the new scattering site, the photon's new radius f' and new comoving frequency v' are 

r' = {r^ + P + 2ri^)^/^ , (26) 

P' = P + 1 (27) 

With the substitutions r' ^ f and u' -^ P, we now repeat Step 2 to find the positions and 
frequencies of the photon at each of its successive scatterings. This loop is repeated until escape 
occurs. 

Step 3: Escape ends the random walk for this photon. We characterize an escaped photon by its 
scaled impact parameter and its "observed" scaled frequency, relative to the source center, 



p = f^l-fi^ (28) 

i^obs = i> -rfi. (29) 



Here the frequency is still defined at the source rest frame and should be divided by (1 + Zs) for 
conversion to the observer's frame. 

Steps 1, 2 and 3 are repeated for as many photons as are necessary to get good statistical 
estimates for the physical quantities of interest. 

In order to relate the results of the Monte Carlo simulations to observable quantities, 
histograms were constructed. Introducing discrete sets of impact parameters pi and "observed" 
frequencies Uj, the final values for the escaping photons can be binned into a two-dimensional 
histogram. The "observed" intensity field I{pi,Uj) is then estimated as proportional to 

^ ("?nl 

where N^j is the number of photons falling into the i,j bin of widths Apj, APj. Similar estimates 
can be made for the intensity /(i5), integrated over frequencies, and /(t'), integrated over the area 
of the plane of the sky. Again, we define these intensities at the source frame. For conversion to the 
observer's frame, /(pj, Vj) should be divided by (1 + Zg)^, and I{p) should be divided by (1 + z^Y 
(since the photon phase space density I{p,u)/h''^ is conserved during the Hubble expansion). 



Although the quantities in equations (^) and ( |27l ) are not directly observed, it is of theoretical 
interest to bin them also into a histogram in order to estimate the value of mean intensity at 
interior points, based on the fact that the number of scatterings per unit volume is proportional 
to the mean intensity. 

A Monte Carlo run using 10*^ photons seemed to give satisfactory results for most of the 
histogram bins. The exceptions are those at small values of the impact parameter or radius, where 
the corresponding bin areas or volumes become quite small, making the statistical errors large. 

The Monte Carlo results are given in Figures 1-5. Figure 1 shows the intensity in the Lya line 
integrated over frequency, 

roo 

i{p)= / i{p,D)dv, (31) 

Jo 

as a function of impact parameter. The intensity has a fairly compact central core; the 

characteristic impact parameter at which the intensity has fallen to half its central value is only 

about O.lr^, an order of magnitude less than the characteristic length r^. Because the intensity 

scales inversely as the square of the characteristic impact parameter, this means that the central 

intensities are two orders of magnitude larger than a simple estimate might have indicated. The 

half-light radius pi/2 (i.e., the radius of a circular aperture containing half of the Lya emission) is 

further out, at 0.63r^, but at that impact parameter the intensity is already down by over an order 

of magnitude from its central value. Since the sky brightness limits the observational sensitivity, it 

may be difficult to infer pi/2 observationally. However, since the scattered light is expected to be 

highly polarized (Rybicki & Loeb 1999), its contrast relative to an unpolarized background could 

be enhanced by using a polarization filter. 
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The angular radius on the sky corresponding to the physical radius of p = 0.1r+ over which 
the scattered Lya halo maintains a roughly uniform surface brightness, is given by 



p_ 



(32) 



where d^ is the angular diameter distance. In an 0,^ = cosmology, ^a = 
2cHq^[Q.mZs + {^M - 2){^/l + nMZs - l)]/[l^if(l + Zsf] (e.g., Padmanbhan 1993), while 
in an 0,\ ^ cosmology the expression is more involved (Edwards 1972, Eisenstein 1997; see also 
Fig. 13.5 in Peebles 1993). We then find for the example of ^m = 0.4, f^A = 0, /iq = 0.65, and 
Zs = 10, that 9 = 15" X {p/70 kpc), with a slightly larger value for a flat J^a = 0.6 cosmology. 
Thus, the Lya luminosity of a typical source at Zg ~ 10 is expected to spread over a characteristic 
angular radius of ~ 15" on the sky. 

Figure 2 shows the total photon emission rate (luminosity) per unit frequency, 

roo 

L{u) = 87r2 / I{p,i))pdp, (33) 



which can be used to get the observed spectral flux of photons F{v') (in photons cm ^ s ^ Hz ^) 
from the entire Lya halo, 

where u = ^^^/{l + z^), and d\^ = dx{^ + -Zs)^ is the standard luminosity distance to the source. 
This flux is expressed in the frame of a local observer. The extra factor of (1 + Zs)^ is due to the 
fact that we evaluate the photon number flux per unit frequency rather than the energy flux - 
which is used in the usual definition of di^. As mentioned before, the observed intensity is 

'^^''^^ = (rT^S' ^''^ 

and similarly the observed integral of the intensity over frequency, 



I{p) iVg 

[l + Zs)^ rl 



m = 77^^^- (36) 



The peak of the function L{D) is seen to be redshifted by roughly 0.44 (recall that we define 
i> to be the negative deviation from line center). The full width at half maximum is AP = 1.75. 
A linear plot of L{i') is presented in Figure 3. Because the shifts are small, the abscissa can be 
equally interpreted as proportional to wavelength shift from line center, so Figure 3 also shows 
the observed lineshape as a function of wavelength. One notes the very extended red wing of this 
line, which goes approximately as l/v^ . Because of this behavior, the centroid of the line is not 
well defined. 

Figure 4 shows the observed intensities as a function of frequency for five different impact 
parameters (logp = 0.0, —0.5, —1.0, —1.5, and —2.0). The curves for \ogp = —1.5 and —2.0 have 
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been smoothed, but still show some fluctuations due to limited Monte Carlo sampling, especially 
at the smallest frequencies. All curves share roughly the same lineshape, and (apart from vertical 
scaling) differ only in the redshift of the peak. The peak redshift is at a minimum, v = 0.19, 
for the central rays, and increases to 0.33 for p = 0.1 and to 1.4 for p = 1. Using appropriate 
integration, these curves could be used to compute the expected lineshape from observations in 
specific circular apertures or slits (e.g., Dey et al. 1998). 

Figure 5 is mostly of theoretical interest and shows the internal mean intensity as a function 
of radius for a variety of frequencies. At the smallest frequencies the curves are reasonably well 



approximated by the diffusion solution (^1|), shown as the dotted curves. As z> increases into 
the free-streaming regime (t' ^ 1), the curves deviate from the diffusion solution, and eventually 
become almost shock-like, with photons piling up near the "causality" surface f = u, consistently 
with the "causal" diffusion equation (18). 

One caveat concerning Figure 5 should be noted. Many of these curves exhibited considerable 
statistical error at small values of f, especially for the largest values of u. In these uncertain 
regions the plots were completed by assuming that all curves approach a constant value at small 
f, determined by the solution at larger f. This extrapolation procedure is certainly justified in the 
diffusion limit, but its validity in the general case needs to be investigated further. 



4. Observational Considerations 

In this section we consider the detectability of the Lya halos. Many of the basic properties of 
the early sources are unknown, and so our estimates are meant for illustrative purposes only. 

The cosmological parameters affect our estimates, but not greatly. In this section, we assume 
f^A = 0, 0,M = 0.4, fifo = 0.05, and /iq = 0.65. Of greater uncertainty is the Lya photon luminosity, 
Na, of high-redshift sources. For the sake of definiteness, we adopt the measured luminosity 
for a known Lya galaxy at Zg = 5.34, which was discovered by Dey et al. (1998). To estimate 
the integral under the observed spectral line profile in their Figure 3, we assume an intensity 
amplitude of ~ 5 /iJy and a rest-frame line width of ~ 2 Ain our cosmological model, this source 
has a luminosity distance of di, = 1.64 x 10^^ cm, yielding a value of 

iV^ = 6.4x 10^^ s-i. (37) 

Let us now assume that a candidate high-reshifted source is identified and its reshift 
determined by some other unscattered lines, such as H^. As a specific example, suppose the source 
is at Zs = 10, so that its observed Lya line falls in the infrared at 1.34/um. From equations (pi) and 
@ we find, 

v^ = 9.8 X 10^^ Hz, 



2.3 X lO^'' cm = 745 kpc. (38) 
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The angular diameter distance is found to be c^a = 3.0 x 10^'^ cm, implying a typical angular size 
of the halo of 0.1 r^/dA = 7.7 x 10"^ = 16". 

The observed integrated brightness near the center of the halo is given by equation (|36|). 
Figure 1 shows that /(O) = 0.2, and so 

4alo(0) = 18 photons cm-2 g-^ sr-i. (39) 

We suppose that a narrow filter is used to observe only in a narrow wavelength band 
containing the Lya line, in order to reduce the sky background as much as possible. From Figure 3 
we infer the Full Width at Half Maximum of the line AzJ ~ 1, implying an observing filter width of 
50 A for a source at z^ = 10. However, one could use narrower filters to take advantage of the fact 
that the monochromatic intensity depends jointly on frequency and impact parameter in a known 
way, as shown in Figure 4. If one wanted to concentrate only on the innermost core of emission 
(out to p = 0.1), then a width Az> = 0.5 would be sufficient, implying a filter width of only 25 
A. We adopt this narrower filter width in our estimate of the sky background. In frequency this 
corresponds to Az^ = 4.2 x 10^^. 

The brightness of the halo is to be compared to the sky background at 1.34 ^m, which is 
obviously minimized in observations from space. Using sufficiently high resolution to eliminate 
contaminating point sources, the sky brightness from space is dominated by interplanetary dust 
(IPD) emission. The local IPD contamination might be reduced by having a spacecraft well above 
the ecliptic plane, but for the moment we shall proceed conservatively and assume that the IPD 
is the primary sky background. Hauser et al. (1998) quote a value for the IPD brightness at 1.25 
^m of ~ 375 nW m~^ sr~^. This implies a sky brightness within our ffiter band of 

/gky = 4.7 X 10^ photons cm~2 s"^ sr"^ (40) 

Thus, the halo brightness is only ~ 4 x 10~^ of the sky background. Given this low value, 
Lya halos would be difficult to observe. However, there are a number of reasons for reasonable 
optimism in the long run for their detection. 



The choice of Lya luminosity in equation (|37|), based on Dey et al. (1998), was made purely 
for definiteness. The observed radiation intensity scales linearly with the value of N^^ and so our 
estimates may be easily modified for any other value. If the values of N^ for some high redshift 
sources were much larger (e.g. due to higher star formation rates or lower dust content), then the 
difficulty in observing the Lya halos would be eased considerably. In this regard we note that even 
if the value given in equation (^) is typical for high-redshift galaxies, it is likely that there will 
be a power-law distribution of observed fluxes (see Figure 2 in Haiman & Loeb 1998b), and some 
sources could have values of N^ larger by factors of ten or more. 

In addition, the interplanetary dust emission could be much lower if the orbit of a spacecraft 
such as the Next Generation Space Telescope (NGST) is chosen so that the telescope spends some 
of its time outside the orbit of the Earth or well above the ecliptic plane. The limiting noise level 
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in this case is provided by the cosmic infrared background. This background is poorly known, but 
Hauser et al. (1998) give an upper hmit at 1.25 fim of i^Iq ~ 75 nW ni~^ sr~^, a factor of five 
less than the value adopted in equation (^0|). Given this sky background, we can calculate the 
minimum exposure time necessary for the detection of the Lya halo signal. Ignoring instrumental 
noise and adopting a detector quantum efficiency close to unity (see http: 
augusta.stsci.edu/ngst-etc/ngstetc.html, for more realistic estimates), we can find the best 
signal-to-noise ratio, S/N, that is attainable after an exposure time t on an NGST telescope of 8 
meter diameter, 

^° ] f^V" (41) 

6 X 1054 s-iy VlO hours/ ^ ^ 

Hence, the halo of a source at Zg ^ 10 which is an order of magnitude brighter in Lya than the 
Dey et al. (1998) galaxy might be detectable. 

Another circumstance favoring such observations is that the Lya halos are highly polarized 
(see Rybicki & Loeb 1999). This polarization is highest at the outermost radii, but is ~ 14% even 
at the core radius of O.lr^. Differencing two maps in orthogonal linear polarization is potentially 
a way to improve the signal-to-noise ratio. More importantly perhaps is that the presence of 
polarization would be a clear signal that the measured halo was a Lya halo of the type described 
here rather than some other, possibly instrumental, effect. 



5. Conclusions 

We have shown that Lya sources before the reionization redshift should be surrounded by an 
intergalactic halo of scattered Lya photons. These sources are expected to appear more spatially 
extended in the Lya line than they are in the continuum to the red of the line. Our numerical 
solution implies that the Lya halo has a roughly uniform surface brightness out to an impact 
parameter p '^ O.lr^, ~ 70 kpc (see Fig. 1). At this impact parameter, the line is broadened and 
redshifted by of order ~ 10^ km s~^ relative to the source (Fig. 3). These substantial broadening 
and redshift signatures cannot be easily caused by galactic kinematics and hence signal the 
intergalactic origin of the scattered line. 

The detection of intergalactic Lya halos with the above characteristics around sources down 
to a limiting redshift at which the neutral IGM ceases to exist, can be used as a direct method 
for inferring the redshift of reionization. Alternative methods are more ambiguous, as they 
rely on the detection the Gunn-Peterson damping wing which might be confused with damped 
Lya absorbers along the line-of-sight (Miralda-Escude 1998; see also Haiman &: Loeb 1999 for 
Zs ^ 7), the detection of the weak damping factor of small-scale microwave anisotropics which is 
an integral quantity depending also on other cosmological parameters (e.g., Hu & White 1997; 
Haiman & Loeb 1998a), or the detection of very faint (and somewhat uncertain) spectral features 
in the cosmic background which is highly challenging technologically (Haiman, Rees, & Loeb 1997; 
Gnedin & Ostriker 1997; Baltz, Gnedin, & Silk 1998; Shaver et al. 1999). 
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Our calculation assumed the simplest configuration of a uniform, neutral, IGM with a 
pure Hubble flow around a steady Lya source. In popular Cold Dark Matter cosmologies, 
the characteristic nonlinear mass scale of collapsed objects at z^ ;^ 10 is ^ 10^ Mq (e.g., 
Haiman Sz Loeb 1998a). A galaxy of total mass 10^ Mg Mq is assembled from a radius of 
AAlMs/^Aih^y'^Kl + Zs)/10]^"'^ kpc in the IGM, which is more than an order of magnitude 
smaller than our inferred Lya halo radius. Similarly, the Hubble velocity at the Lya halo radius 
is larger by more than an order of magnitude than the characteristic velocity scale of nonlinear 
objects at these redshifts. Hence, our simplifying assumptions of a smooth IGM immersed in a 
Hubble flow are likely to be satisfied on the Lya halo scale. Modest corrections due to the density 
enhancements and peculiar velocities in the infall regions around sources might, however, be 
necessary. More importantly, the ionization effect of a bright quasar on its surrounding IGM might 
extend out to the scale of interest. In such a case, the intensity distribution of the Lya halo will 
depend on the spectrum and luminosity history of the ionizing radiation emitted by the source, 
which determine the neutral fraction as a function of radius around it. This "proximity effect" 
might be important for quasars but less so for galaxies whose ultraviolet emission is typically 
strongly suppressed beyond the Lyman limit due to absorption in stellar atmospheres and in the 
interstellar medium. Other changes to the halo intensity and polarization profiles might result 
from short-term variability (on ^ 10^ years) or anisotropic Lya emission by the source. These 
complications could be easily incorporated into our Monte Carlo approach, for particular source 
parameters. 

Detection of the predicted Lya halo might become feasible over the next decade either 
with larger ground-based telescopes or with the Next Generation Space Telescope^ (NGST). For 
Zs ~ 10, the entire Lya luminosity of a source is typically scattered over a characteristic angular 
radius of ~ 15". The Lya halo is therefore sufficiently extended to be resolved along with its 
tangential polarization, as long as its brightness exceeds the fluctuation noise of the infrared 
background. 

Recently, five sources have been photometrically identified to have possible redshifts of z ;^ 10 
in the Hubble Deep Field South observed by NICMOS (Chen et al. 1998). We emphasize that 
even just a narrow-band photometric detection of the scattered Lya halos (see Fig. 1) around 
sources at different redshifts would provide invaluable information about the neutral IGM before 
and during the reionization epoch. On sufficiently large scales where the Hubble flow is smooth 
and the IGM is neutral, the Lya brightness distribution can also be used to determine the values 
of the cosmological mass densities of baryons and matter through equations (^), ^ and the 
angular diameter distance relation. Thus, in addition to studying the development of reionization, 
such observations could constrain fundamental cosmological parameters in a redshift interval that 
was never probed before. 



^NGST is the successor to the Hubble Space Telescope which is planned for launch over the next decade. For 
more details, see |ittp://ngst. gsfc.nasa.gov/ . 
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Scattering of resonance line radiation gives rise to polarization at a level depending on the 
atomic physics of the transition and on the geometry of the scattering events. Using an extension 
of the Monte Carlo technique used here, Rybicki & Loeb (1999) showed that most photons 
with the largest impact parameters suffer a nearly right-angle last scattering, and the degree of 
tangential polarization is large for such photons, of order tens of percent (asymptotically 60%). 
Measurement of this polarization could be used as an independent check on the model parameters, 
and might even provide a powerful way of identifying early objects. 

We thank Chris Kochanek and Chuck Steidel for useful discussions. This work was supported 
in part by the NASA grants NAG5-7768 and NAG5-7039 (for AL). 
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Fig. 1. — Integrated intensity, I{p) = /q°° IdD, versus impact parameter, p. The projected half Hght 
radius is denoted by Pi/2- 
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Fig. 2. — Monochromatic photon luminosity, L{D) = Svr^ /q°° / pdp, as a function of frequency 
redshift, P. 
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Fig. 3. — Same as Figure 2, but using linear scales. 
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Fig. 4. — Lya line profiles. Curves of intensity, I{p, D), versus frequency shift to the red of resonance, 
P, at different values of the impact parameter, p. 
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Fig. 5. — Curves of the monochromatic mean intensity, J(r, i>), versus radius, f, at different values 



of the frequency redshift, u. Note that in agreement with the "causal" diffusion equation (18), 
^ ^ /^ ^^) where / = 1/3 in the diffusion regime (f <^ 1) and / = 1 in the free-streaming regime 
(f 3> 1). The dotted lines show the analytic solution (ETl) for the diffusion regime. 



